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Abstract 

We construct the potentials that describe the spectrum and decay of elec- 
tromagnetic bound states of hadrons, and are consistent with ChPT. These 
potentials satisfy the matching condition which enables one to express the 
parameters of the potential through the threshold scattering amplitudes cal- 
culated in ChPT. We further analyze the ambiguity in the choice of the short- 
range hadronic potentials, which satisfy this matching condition. 
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I. INTRODUCTION 



The energy spectrum and decays of electromagnetic bound states of strongly interacting 
particles - so called hadronic atoms - have recently been measured by several experimental 
collaborations. These measurements yield an extremely valuable piece of information about 
the behavior of QCD at a very low energy, which is hardly accessible with a different ex- 
perimental technique. In particular, the measurement of the tt+tt" atom decay width by 
DIRAC collaboration at CERN |I| which will result in the determination of the difference 
Oq — 0,2 of the S'-wave vrvr scattering lengths at a 5 % precision, would allow one to directly 
test the large/small condensate scenario of chiral symmetry breaking in QCD with two fla- 
vors. Further, the Pionic Hydrogen collaboration at PSI intends to extract the S'-wave t^N 
scattering lengths from the ongoing measurement of the spectrum and transition energies 
between 3p — Is levels in pionic hydrogen at a 1 % accuracy which is unique in hadron 
physics 0]. This will yield a more precise value of the ttNN coupling constant and of the 
T^N o"-term. Finally, the DEAR collaboration 0] at the DA$NE facility plans to measure 
the energy level shift and lifetime of the Is state in K~p and K~d atoms - with considerably 
higher precision than in the previous experiment carried out at KEK for K~p atoms. 
It is expected that this will result in an accurate determination of the J = 0, 1 S'-wave 
scattering lengths. It will be a challenge for theorists to extract from this new information 
on the KN amplitude at threshold a more precise value of e.g. the KN cx-term and of the 
strangeness content of the nucleon. 

In order to fully exploit the high-precision experimental data, it is imperative to design 
the theoretical framework for the analysis of these data which would describe this sort of 
bound states in the accuracy that matches the experimental precision. In practice, this 
would mean that the Deser-type formulae [§] which are used to extract the strong scatter- 
ing lengths from the measured values of the energy shift and width of hadronic atoms, do 
not provide the necessary precision and should be corrected to accommodate strong and 
electromagnetic isospin-breaking corrections which typically amount up to a few percent. 
Since these corrections arise in a result of a complicated interplay of strong and electromag- 
netic effects in the bound-state observables, some dynamical input on strong interactions is 
needed for their evaluation. Starting from the problem of the calculation of hadronic 
atom observables has been analyzed within the potential scattering theory approach (see, 
e.g. [[f|-[TO[|), where the strong interactions are described by the short-range "hadronic" poten- 
tials. In order to possess a predictive power, certain assumptions about the strong potential 
should be made within this approach - it is usually assumed that the short-range strong 
potential does not violate the isospin symmetry and describes purely hadronic data at low 
energy. Isospin-breaking corrections are then calculated by taking into account the kinemat- 
ical effect due to the mass difference of charged and neutral particles, including Coulomb 
interaction, vacuum polarization and some higher-order pure QED effects. 

The issue of "purification" of the experimental data with respect to strong and electro- 
magnetic isospin-breaking corrections arises in the context of the the low-energy hadronic 
scattering processes as well. This problem is closely related to the problem of the hadronic 
bound states, which were discussed above. Two distinct approaches are available in the liter- 
ature to address the issue in the scattering sector: the dispersion relations approach |TT| , p!2 



and the potential scattering theory approach (see, e.g. [p!^-[l5|). The results obtained within 
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these two approaches are used e.g. in the recent analysis of the low-energy vrA^ scattering 
data Hill. 



Recently, the problem of isospin-breaking corrections in hadronic atom observables has 
been addressed in the framework of ChPT |[17|^2| (see Ref. for the comparison of 



different approaches). In a result of these investigations, a systematic discrepancy between 
the results obtained within the potential scattering theory approach and the field-theoretical 
approach, has been discovered for the case of the 7r"^7r~ atom decay width and the 7r~p atom 
ground-state energy. In particular, the prediction for the vr+vr" atom decay width |^-|^ 
depends on the potential used, and on the choice of the free Green function in the Schrodinger 
equation. Earlier calculations |^ have predicted the isospin-breaking corrections with the 
opposite sign and with the same order of magnitude as in ChPT. In the later paper 
including some kinematic relativistic corrections in the free Green function, the authors 
obtain a result for the lifetime r = 3.0 x 10~^^ sec which is still not compatible with 
the existing prediction of ChPT: for a fixed value of 09 — 02, the lifetime equals to r = 
2.9 X 10^^^ sec with an accuracy of 1.2 % [Q. The latest prediction of the potential model for 
the lifetime is numerically close to that of ChPT 0. However, no compelling physical reasons 
are available in favor of the conclusion, that the latest version includes all relevant isospin- 
breaking corrections from the underlying theory. The potential-approach based prediction 
for the energy shift of the 7i~p atom [|10[ differs by a factor 2 from the order-of magnitude 
estimate carried out in ChPT at O(p^) [^, and quotes the value 0.5% for the systematic 



error which is 4 times smaller than the corresponding value obtained from ChPT. Even if 
the values of the low-energy constants (LECs) in ChPT turn out to be such (see Ref. p3| for 
the evaluation of these constants in the quark model), that the prediction for the correction 
in ChPT comes numerically close to that of the potential approach, still the systematic error 
is underestimated in the potential approach. Of course, the large systematic error in the 
prediction of ChPT comes from the LECs which, albeit poorly known in the vrA^ case, are 
still present and can not be simply disregarded. In addition, we would like to stress here 
that there is no physical reason why the LECs in ChPT must have the particular values 
which reproduce the results of the potential approach. 

The discrepancy mentioned above does not come to a surprise: the potential which is used 
to calculate the bound-state observables, does not include a full content of isospin-breaking 
effects in QCD+QED. It has been demonstrated (see, e.g. ||2^) that the isospin-symmetric 
short-range potential which has been used so far, can not fully accommodate the effect of 
the direct quark-photon interaction, as well as the isospin-breaking effect due to the explicit 
dependence of the scattering amplitude on the quark masses which is governed by the chiral 
symmetry. These effects turn out to be dominant in the isospin-breaking part of vrvr and vrA^ 
interactions at threshold, that leads to the above-mentioned discrepancy. Moreover, since 
the potential approach makes use of the same ideology in the calculation of the isospin- 
breaking corrections to hadronic scattering amplitudes in the low-energy region, the same 
criticism is applicable in this CclSG clS well - in fact, a significant discrepancy between the 
results of the potential approach and those of the ChPT in the analysis of tcN scattering 
process has already been reported p5| . 

The aim of the present paper is to explicitly demonstrate, how the potential model can 
be made compatible with ChPT - in what concerns the calculation of the isospin-breaking 
corrections to the observables of hadronic atoms. This construction will help one to finally 
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resolve the long-standing puzzle and to determine the status of the potential approach-based 
calculations in the analysis of the hadronic atom data. Moreover, in the view of the fact that 
the potential approach is widely used to evaluate the isospin-breaking effects in scattering 
processes, it is important to note that the same construction allows one to constrain the 
potential from comparing the threshold amplitudes in the field theory and in the potential 
approach. In this vein, one may hope to reduce the ambiguity in the prediction of the isospin- 
breaking corrections above threshold. The ultimate goal of our investigation is to derive the 
short-range hadronic potentials - including the part which violates the isospin symmetry - 
from ChPT. Although extensive investigations address the similar issue in the context of 



purely strong potential (see, e.g. p6|-p8|), to the best of our knowledge, the comprehensive 



study of isospin-breaking effects is not available in the literature. 

The layout of the present paper is follows. In Section ^ we briefly review the derivation of 
the potential in the absence of isospin breaking. Section |ITI| deals with the simplified case of 
one-channel scattering, whereas in Section ^ we address the inclusion of the isospin-breaking 
effects in full generality (multichannel scattering, relativistic corrections). In Section we 
compare the results obtained within ChPT and existing potential models. Finally, Section |V| 
contains our conclusions. 



II. STRONG POTENTIAL 

For demonstrative purposes, below we shall briefly review the basic idea which lies in 
the foundation of the present approach. Assume that one starts from a given relativistic 
field theory, and is aimed at the construction of the potential which, used in the Lippmann- 
Schwinger equation, produces the same ^-matrix elements at low energies, as the initial 
relativistic field theory. In order to achieve this goal, it has proven very useful, at an inter- 
mediate stage, to construct the non-relativistic effective field theory that correctly describes 
the low-energy behavior of the initial relativistic theory. To ease the consideration, we treat 
here the simplest case of scalar self-interacting field. Electromagnetic effects are assumed to 
be absent. In the low-energy domain (at energies much less than the mass of the particle), 
the physics can be described on the basis of the non-relativistic Lagrangian 

(2.1) 

where M is the physical mass of the particle, u /\ v = u{Av) + {Au)v, and 0(a;) is the 
non-relativistic field operator 

0(0, x) = I duip) e^^- a(p) , [a(p), at(k)] = (27r)3<53(p _ k) , rfz/(p) = ^ . (2.2) 

The Lagrangian (|2.1| ) contains an infinite string of local operators with an increasing mass 
dimension, which describe the scattering process (jxp 4>4>- The operators which contain 
powers of space derivatives are suppressed by the corresponding power of M which is the 
only heavy scale available. Thus, the contribution of these operators becomes increasingly 
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suppressed^ at a small momenta |p| << M. Note that higher orders in time derivative dt 
can be systematically eliminated by the use of the equations of motion, without altering the 
S*- matrix elements (see, e.g. p9|). 



The detailed discussion of the perturbation theory based on the Lagrangian (|2.1|) can be 



found in Refs. p2| , p9| - below, we only provide the sketch of the procedure. The S'-matrix 
element for the scattering process in the non-relativistic theory is given by a sum of tree 
diagrams, plus any number of s-channel bubbles, plus mass insertions in the external and 
internal lines which come from the higher-order terms in the kinetic part of the Lagrangian 
(see Fig. |T|). The free propagator of the scalar field is given by 

'»)^''°)l°)-/|^ M + pV(2M)-pO-. - P-^' 

where (j) denotes the free field. To ease notation, we omit —ie term in the following. The 
elementary building block in the calculation of the diagram with any number of strong 
bubbles in the CM frame is given by 

j( pO) = f I I = !^ (M( - 2M))^/^ 

^ ' 7 (27r)^i M + 1V(2M) -PO + /o M + IV(2M) -/o 47r ^ ^ " 

(2.4) 

at — s> 4, P° > 2M. The scattering matrix elements are obtained by putting P° = 2w(p) 
where w(p) = (M^ + p^)^^^ is the energy of the particle in the CM frame. At threshold, 
in dimensional regularization, the integral is purely imaginary and is proportional to i|p| 
where p is the three-momentum of the particle in the CM frame. The integrals containing 
derivative vertices and/or mass insertions can be considered analogously. According to the 
standard power-counting, their contribution will be suppressed by the corresponding power 
of |p|. 

The matching condition for the relativistic and non-relativistic scattering am- 
plitudes in the CM frame reads 

T^(p, q) = Aw\v) T^''(P, q) , |p| = |q| • (2.5) 

This matching condition is understood as follows. At threshold, both relativistic and non- 
relativistic amplitudes can be expanded in powers of CM momenta p and q 

T' = fo + |p|/i' + pV2 + Pq fi + O(p^) , I = R,NR. (2.6) 

The matching condition then gives the relation between the coefficients of the expansion 
in the left and right-hand sides of Eq. ( p.5| ). In the lowest order in p, the non-relativistic 
scattering amplitude is completely determined by the tree diagram containing the coupling 



^We consider only the case of the so-called "natural EFT" where such an estimate for the size of 
the couplings is justified, and one may straightforwardly apply the perturbation theory to calculate 
the ^-matrix elements. 
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qq (Fig. |T]a) - bubbles, derivative vertices and mass insertions all give the contributions that 
vanish at threshold. Consequently, the matching condition in the lowest order in p reads 



{2\ fgo = r^(0,0) 



(2.7) 



where in the right-hand side the S'-wave scattering length in the relativistic theory appears 
(here 2! is the symmetry factor for the identical particles in the initial and final states). 
Note that in the matching, we have not used the perturbative expansion in the relativistic 
theory, so the above relation is valid "in any order in strong coupling constant". Further, 
using the same technique, the coupling gi can be related to the S'-wave effective radius and 
the S'-wave scattering length, the coupling g2 can be related to the P-wave scattering length, 
and so on. 

In the following discussion, we shall first neglect derivative couplings, as well as mass 
insertions p^/{8M^) ■ ■ ■. These will be considered in Section |V[ The non-relativist ic scat- 
tering amplitude in this case is given by 



T 



NR, 



p,qj 



4^0 



1 



27r \t>\ 



(2.8) 



Now it is straightforward to check that all strong bubbles with non-derivative vertex g^ can 
be resummed with the help of the following Lippmann-Schwinger equation 



r(£; P, q) = V(p, q) + ^ / (gj 1'(P. 1) nE; 1, q) 



(2.9) 



and 



-T 



NR, 



|q| = |p|, £;=q2/M 



|q|=|p| 



(2.10) 



where V^(p, q) = —{21)'^ go, and we have used dimensional regularization to regulate UV 
divergences in physical dimensions (i — > 3. Note that the operator T{E) acts in the Hilbert 
space of vectors |p), whereas T^^(p, q) are the matrix elements of the scattering operator 
j^NR ^Yie Fock space between two-particle states |p, — p) = a^(p)a^(— p)|0), with the 
CM momentum removed (for details, see |jT9|j2^). Below, to ease discussion, we shall not 
explicitly distinguish between the state vectors in the two spaces. 

Now we address our main task of deriving the potential from initial relativistic field 
theory. In principle, Eq. ( |2.9| ) already solves this problem - in the dimensional regularization. 
The potential has then a 5-type singularity in the coordinate space 



V(r,r') 



^ipr— iqr 



{27rY (27r) 



V{p,q) = -i2\fgoS\r 



r') 6''(r) 



(2.11) 



In the conventional scattering theory it is customary to deal, however, not with the pointlike 
interactions in the dimensional regularization, but rather with the smooth potentials in the 
position space in the physical space dimensions d = 3. We may accommodate this feature by 
smoothing the potential which is obtained from field theory, and simultaneously adjusting 
the coupling g^ so that the scattering amplitude still obeys the matching condition (^.71). 
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We demonstrate this procedure for the simplest case when one assumes the separable ansatz 
for the resulting potential 



r(r,r')^V0(r,r') 



-(2!)2^7o(/?)^(r)t;(r') 



f r 



-/3r 



47rr 



/52 + p2 



(2.12) 



where (7o(/?) stands for the potential strength, and j3 is the potential range parameter (we 
anticipate that, in a result of matching, §^{[3) will depend on {3). The scattering amplitude 
in the case of the separable interaction is given by 



T(E;p,q) = -(2!)2^o(/3)t^(p)t^(q) 



l + 2!^7o(/5) 



(27r)3 E-F/M 



-1 



The matching condition for the scattering amplitude at threshold then gives 



AM' (2!f^?o(/?) = 



r^(o,o) 

T«(0, 0) 



1 + ^ 

^ 647rM 



(2.13) 



(2.14) 



Of course the choice of the separable ansatz is not unique - one may, e.g. use the local ansatz 
V(r, r') — > 5^{y — v') Va(r) where j3 denotes a range parameter for the potential. One may 
further choose any form for the potential Va(r), e.g. use the double square well potential [0. 
In the absence of the analytic solution in this case one will have to solve the matching 
condition for the coupling constant numerically. The generalization to higher-order terms 
in expansion, as well as for higher partial waves is straightforward. 

To summarize, we note that, according to the point of view adopted through the present 
paper, the smooth potentials in the position space are nothing than a particular regulariza- 
tion of the genuine pointlike interactions between particles, which stem from the underlying 
field theory. The shape of the potential thus bears no important physical information. The 
physical input from the underlying theory is contained in the couplings, and enters through 
the matching condition, which states that the fixed number of terms in the effective-range 
expansion of the scattering matrix elements coincide in the potential approach and in the 
underlying relativistic field theory. This is enough for both approaches to describe the same 
physics in the vicinity of the scattering threshold. 



III. ISOSPIN-BREAKING EFFECTS IN THE ONE-CHANNEL CASE 

In this section, we discuss the inclusion of the electromagnetic effects in the scheme 
considered above. Once the photons are included, a qualitatively new feature in the theory 
emerges: particles with opposite charges can be bound by Coulomb force at a distances 
much larger than a typical range of strong interaction. The observable characteristics of 
this sort of bound states - energy levels and decay widths - receive contributions from 
strong interactions. This fact enables one to extract strong scattering lengths from precise 
measurements of these characteristics, if a systematic quantitative description of such a 
bound states is provided on the basis of the underlying theory (see, e.g. |T9|-p2|). 

In order to extract the parameters of the strong interactions from the experimental quan- 
tity which contains contributions from both - electromagnetic and strong - interactions, one 
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has to say explicitly, how such a splitting can be systematically performed in the framework 
of the field theory. As a simplest example, one should be able to disentangle electromagnetic 
and strong contributions in the mass of a particle which takes part in both interactions. In 
general, the issue has proven to be rather subtle - the splitting is convent ion- dependent. 
Further we do not discuss this question, since it forms a separate subject for the investi- 
gation. Note that in all examples considered below, one may choose the convention in the 
underlying relativistic theory so that the mass of the particle in the "purely strong" world 
coincides with the mass of charged particle in the real world, and the explicit prescription 
for the splitting of the scattering amplitude into the strong part and the electromagnetic 
correction can be provided. Below, we shall follow these conventions. 

The question which we investigate here can be now formulated as follows. We study the 
observables of electromagnetic bound states - energy and decay width in a theory where both 
strong and electromagnetic interactions are present. In the leading order in fine structure 
constant a, the relation between these observables and the strong scattering lengths is 
universal 0. Nontrivial interplay between electromagnetic and strong effects occurs at the 
first non- leading order in a. Our aim is to evaluate these observables at the first non- 
leading order, and to construct the short-range hadronic potential which, when used in the 
Schrodinger equation, leads to the same values of these observables at the same order in a. 



A. Energy-level shift in field theory 

The relativistic field-theoretical Lagrangian of the model which is considered in the 
present Section, describes the charged scalar particle with the mass M, interacting with 
photons. In addition, the Lagrangian includes arbitrary self-interactions of the scalar par- 
ticle. In the theory, there are two-particle loose bound states with zero total charge, which 
are formed mainly by the static Coulomb interaction. In the lowest order in a, the energy 
of the ground state coincides with the ground-state energy in the pure Coulomb potential: 
Eq = 2M — i a^M, and, as was mentioned above, there are higher-order corrections to this 
value, caused both by the electromagnetic and strong interactions. Our goal is to calculate 
the energy of the ground state up to and including 0{a^), where both these sources con- 
tribute. It is both convenient and conventional to split the ground-state energy at 0{a'^) in 
the electromagnetic part and strong shift, according to pi]| , pT 



The expression for the electromagnetic part can be easily obtained, adapting formulae from 
Ref. to the scalar caseQ 

1 S 1 

= Eo--^Ma^ + ^ M^{r'')a^ + o(«^) , (3.2) 

where (r^) denotes the charge radius of the particle. Assuming (r^) = 0, we reproduce the 
result for the ground-state energy-level shift at O(a^), given in Ref. [EO. 



^It suffices to substitute A = | (r^) in Eq. (4) of Ref. |21] and evaluate energy shift according to 
the formulae (13) and (14) at nip = = M. Vacuum polarization contribution is omitted. 



In the following, we shall concentrate exclusively on the strong shift at 0{a'^) which is 
given by the second term in Eq. (|3.1| ). 



1. Non-relativistic Lagrangian 

The non-relativistic Lagrangian which is sufficient to calculate the strong shift at O(a^), 
consists of the non-relativistic kinetic term, static Coulomb interaction, and the short-range 



interaction described by the local four-particle non-derivative vertex |21 



Co = Yl ^\ i^dt -M + 

i=± \ 



A 



2M, 



+ 



t 



/:5 = c(0^,0L)(</)+0_) + -- - , (3.3) 

where the non-relativistic fields (j)± describe the pair of charged particles. All higher-order 
derivative terms that one may write here, do not contribute to the strong shift at the accuracy 



we are working (we refer the reader to Ref. to more details). Further, the coupling 



constant c in Eq. ( p.3| ) is determined from matching the relativistic and non-relativistic 
scattering amplitudes at threshold (see below). Generally, one may write c = Co+aci+0{a'^), 
where Cq is obtained by matching to the strong relativistic Lagrangian in the absence of 
electromagnetic interactions. 

On the potential scattering theory language, the non-relativistic Lagrangian ( |3.3| ) de- 
scribes one-channel scattering problem for Coulomb + strong interactions, in the channel 
with total charge Q = 0. The only difference is, that now the strong coupling c - through 
matching to the underlying relativistic theory - depends on a as well. This is the way how 
the full content of isospin-breaking effects in the initial theory enter in the non-relativistic 
framework. 



2. Hamiltonian framework and the energy of ground state 

It is convenient to use the Hamiltonian formulation of the non-relativistic theory. The 
Hamiltonian, derived from the Lagrangian ( |3.3| ) with the use of canonical formalism, is given 

by 

H = Ho + Hc + Hs, Hr = 1 c^=^x7^^(0,x), r = 0,C,^ 

Hc = ^7ra (0L0_)A-i(0V0+) 

7^5 = -c(0Wl)(0+0_) (3.4) 
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The scattering states in the non-relativistic theory are |P,p) = a+(pi)aL(p2) |0), where 
P = Pi + P2 and P = I (pi — P2) are the CM and relative momenta, respectively. The CM 
momentum is removed from the matrix elements of any operator R in Fock space by using 
the notation 



(P,q|i?(z)|0,p) = (27r)^5^(P)(q|r(z)|p), (3.5) 

The operator r{z) acts in the Hilbert space of vectors |p), where the scalar product is defined 
as the integral over the relative three-momenta of particle pairs 

(q|p) = (27r)3 53(q_p), 1 = | dz/(p)|p)(p| . (3.6) 



The Hamiltonian ( p.4|) at c = has an infinite number of pure Coulomb discrete eigen- 
states. The ground-state eigenvector in an arbitrary reference frame is given by 

1^0, P) = / ^v^o(q)|P,q), (3.7) 
where ^E'o(p) is the Coulomb wave function in the momentum space 



[p2 + 7')' ' ^~ 2 



(p|M/o) = ^o(p) = fy^^, i = -aM. (3.8) 



The eigenstate |\E'o, P) satisfies the following bound-state equation 

(Ho + He;) |vl/o, P) = (i?o + ^) 1^0, P) , Eo = 2M-^ Ma"" . (3.9) 

The resolvent in pure Coulomb theory is defined by Gc{z) = {z — Hq — Hc)^^. In the 
CM frame this resolvent is given (cf with Eq. (|3.5|) ) by the Schwinger's representation ^1 



(27r)353(q_p) 1 

(q|gc(2;)|p) = 



E-(f/M E-qVM |q-p|2 ^-pVM 

^ Anar]IiE; q, p) ^ , (3.10) 



with 

cc ^ dec 

I{E- q, p) = _ ^^2^ ^ ^2/^2(1 _ 2.)2(^ _ q2/M)(E - pVM)] ' ^^'^^^ 

where r] = \ a {-E/M)-^/^ and E = z - 2M. 

When c 7^ 0, the bound-state poles in the full resolvent G{z) = {z — Hq — He — Hs)^^ are 
shifted from their original Coulomb values. Applying Feshbach's formalism we obtain 
the equation for the shifted ground-state pole position p9|j21|j2^j29| . 

z-Eo-i^o\riz)\^o)=0 (3.12) 

where t{z) obeys the Lippmann-Schwinger equation 
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Tiz) = hs + hsgcizHz) , (q|gc;(^)|p) = (q|gc(^)|p) - ^"^'^^"^^^"'^^ - (3.13) 

Here gc(^) denotes the pole-subtracted Coulomb Green function, and (q|h5|p) = — c. 

The bound-state equation ( |3.12| ), together with Eq. ( |3.13D , can be solved iteratively, 
along the similar lines as in Refs. |[T9|j21|j2^ . We use the dimensional regularization both for 
ultraviolet and infrared divergences. The key observation is, that in dimensional regulariza- 
tion, the subsequent terms in the iterative solution of Eq. ( p.l3|) for the quantity t{z) are 
suppressed by powers of a (roughly, one may count h5 ~ 1, gc ~ a), so that only a finite 
number of iterations survives in a given order in a. At 0(a;^), we obtain 



a 



z-Eo = — cn--— c A(^) + ln— + 21na-3 +»(«)> (3-14) 



where 

A(/.) = if^Y'' - r'(l) - ln47r) , (3.15) 

and (i — > 3 in physical space dimensions. 

We define the renormalized coupling constant c^lfi) as 

c = co + aci + 0{a^) = c'-(/i) + -— [c^'(/^)]'A(/i) + 0{a') . (3.16) 

OTT 

Note that both cq and ci can be assumed to be real. At the lowest order in fine structure 
constant, the imaginary part of c is determined by the process 0+0- —>■ 27 in the relativistic 
theory and starts at O(a^). 

Expressed in terms of c''(/i), the energy-level shift does not contain the ultraviolet diver- 
gence 

AE = Rez-Eo = Rec'^(/i) <^ 1 - Rec7/i) In — + 21na - 3 + o(a^) . 

Stt [ Stt V a* /J 

(3.17) 

It is easy now to check that AE does not depend on fi at 0{a'^). 



3. Matching to the relativistic amplitude 

Above, we have expressed the strong energy-level shift of the ground state in terms of 
the renormalized coupling constant c'"(/i) in the non-relativistic Lagrangian. In order to 
have connection to the initial relativistic theory, one has to express this coupling - through 
the matching of relativistic and non-relativistic scattering amplitudes at threshold - via the 
parameters of the relativistic Lagrangian. From Eq. ( |3.17| ), one may conclude that c''(/i) 
should be known at 0(a), so that it suffices to match the amplitudes calculated at the same 
accuracy. Further, note that, in order to obtain the strong shift, one has to perform matching 
for truncated amplitudes obtained by discarding (both in the relativistic and non-relativistic 
theories) all diagrams that are made disconnected by cutting one photon line: T = T — T^x, 
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where Tex stands for the one-photon exchange contribution pT| . The contribution from Tex 
is contained in the electromagnetic shift (|3.2| ). 

The matching condition for truncated amphtudes in the CM frame is analogous to 
Eq. ( p.5|) in the strong sector 

f«(q,p)=4«;2(p)f^^(q,p). (3.18) 

In the non-relativistic theory, the full scattering amplitude at 0{a) in the vicinity of 
threshold is determined by the diagrams depicted in Fig. |^ which are calculated by using the 
Lagrangian ( p. 31) . The truncated amplitude T^^ is obtained by discarding the contribution 
coming from diagram in Fig. ^ which corresponds to the exchange of the Coulomb photon. 
Further, the non-relativistic amplitude contains the (divergent) Coulomb phase which is 
removed. Finally, the expression for the real part of the non-relativistic truncated amplitude 
at threshold is given by 

B^^ 2lr)l 
_ ^1 , nlyR^„ ^IPI 



Re ^«(q, p)] 

where 



IqNIpHo IpI m 



+ 5rin4^ + ^^" + 0(|p|), (3.19) 



Br = ^ Rec^(/i) + o{a) , = [Rec^(/.)]^ + o{a) 



^^^ = Rec'-(/i)|l-^Rec"(/i) (^In ^ - l) | + o(a) , (3.20) 



and 



The scattering amplitude in the relativistic theory contains the diagrams that can be made 
disconnected by cutting one virtual photon line Fig. |^, and the rest which, by definition, 
coincides with truncated amplitude Fig. ^d. After the removal of Coulomb phase, the trun- 
cated relativistic amplitude, like its non-relativistic counterpart, in the vicinity of threshold 
contains singular terms 



Re [e-2^°^cT^(q, p)] 



' B^ln^+A'' + 0{\p\). (3.22) 



iqi=ipi^o IpI m 

The matching condition ( |3.18| ) for the regular parts of the amplitudes at threshold gives 

A^ = AM^A^^. (3.23) 
Using Eq. ( |3.20| ), one may relate c^{fi) with 

1 ,p r. a ,p / 



Finally, substituting this relation into the expression for the strong energy- level shift, we get 

AE = ^« |i _ ^ ^« (In a - 1) j + o(a^) , (3.25) 

327r I loTT J 

where reference to the non-relativistic theory has completely disappeared - the strong shift 
is expressed in terms of the regular part of the relativistic scattering amplitude at threshold. 
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4- Matching condition for the non-relativistic couplings 



As pointed out above, we implicitly assume the prescription which enables one to split 
the scattering amplitude into the strong piece and the electromagnetic correction which is 
proportional to a 

A'^ = A^ + aA^ + Oia^) . (3.26) 
Then, from the matching condition ( p.24[ ), for the non-relativistic coupling we get 



where c = cq + aci + 0(q;^). Note that the equations for determining the energy of the 
bound state are formally the same in the non-relativistic effective Lagrangian approach and 
in the non-relativistic potential model. This again, as in the purely strong case, allows 
one to interpret the local 4-particle interaction in the Hamiltonian ( p.4|) as the contact 
potential V{r) = —c6'^{r) with the strength determined from matching to the relativistic 
field theory (|3.27|) . If one uses the dimensional regularization throughout, this is a perfectly 



consistent interpretation - we have seen that the ultraviolet divergence contained in c cancels 
with the divergence arising in the bound-state calculations. We see also that the coupling c 
contains the piece proportional to a which is related to the corresponding piece a in the 
expression of the relativistic amplitude. This means that, in general, one can not assume 
that the contact term corresponds to purely strong interactions. 

Below, we shall demonstrate that, in a complete analogy with a purely strong case, one 
may also construct conventional smooth short-range potential that reproduces the answer 
for the energy of the ground state, obtained within the field-theoretical approach. 



B. Energy-level shift in potential model: universality 

Below, we consider the non-relativistic potential model with the potential given by a 
sum of Coulomb and short-range parts 

{q\v\p) = v{q,, p) = + u(q, p) = (q|hc;|p) + ^(q, p) . (3.28) 

|q-pr 

One possible interpretation of this potential is to be the ultraviolet-regulated version 
of the contact potential introduced in the previous section, to which it reduces in the 
local limit M(q, p) — c. According to the discussion above, the short range poten- 
tial, in general, should contain the strong part, as well as the piece proportional to a: 
^(QjP) = ^o(q5P) + Q;Mi(q, p) + O(a^). Our goal is to construct the short-range potential 
u(q, p) which reproduces the result (|3.25|) for the energy shift of the ground state, obtained 
within field theory. One would expect, that to this end it is necessary to perform match- 
ing of scattering amplitudes at threshold according to the same matching condition (|3.18 ) 
as in field theory. This should result into the matching condition between the relativistic 
threshold amplitude A^ and the short-range potential M(q, p) which resembles the matching 
condition (|3.27|) for the non-relativistic couplings. 
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Below, we shall demonstrate that the above conjecture indeed holds for the model consid- 
ered here. To be more specific, one has to prove the universality of the relation ( p. 25 ): The 



relation between the energy shift of the ground state and the threshold scattering amplitude 
at next-to -leading order in a is the same in field theory and in the potential model. 

In order to prove this statement, we look for the ground state pole in the scattering 
matrix which obeys the Lippmann-Schwinger equation 

T{z) =v + v go{z)T{z) , (q|go(.)|p) = ^^^^'f^^^ • (3.29) 

The relation of T{z) to the scattering amplitude T^^(q, p) is given by Eq. ( p.lO|) . The 
shift in the ground state pole position in the potential model is again given by Eq. ( p. 121) , 
where t{z) now satisfies the equation ( |3.13| ) with (qlh^lp) = M(q, p). 

An important remark is in order. In the non-relativistic effective field theory, we have 
used dimensional regularization to handle ultraviolet divergences. In this scheme, one may 
solve the equation for t{z) iteratively, since, as it is easy to see, higher-order iterations 
amount to higher-order contributions in a to the bound-state energy. This is no longer 
the case for the short-range potential M(q, p): iterations corresponding to the free propa- 
gation of particles should be summed up in all orders. In order to do this, note that the 
Schwinger's representation for Coulomb Green ( p.lO| ) function suggests its decomposition 



into 0-Coulomb, 1-Coulomb and many-Coulomb pieces. Consequently, 

gc{z) = g°-^(^) + g'-^'{z) + g^-'^iz) = go{z) + 5g{z) . (3.30) 

Whereas in the above decomposition 6g{z) which contains the exchanges of at least one 
photon can be still treated as a perturbation, one has to sum all iterations of go{z). To this 
end, we introduce the auxiliary scattering matrix 

t{z) =u + ugo{z)t{z) . (3.31) 

The relation between t{z) and t{z) is given by 

T{z)=t{z)+t{z)6g{z)T{z). (3.32) 

The partial-wave expansion of the potential and scattering matrix is given by 

(q|n|p) = 47r ^ lLA/(q) UL{q,p) >^lm(p) > 

LM 

(q|t(2;)|p) =4nJ2 >LA/(q) tL^z; q,p) Y^^,{p) , 

LM 

iq\T{z)\p)=AnJ2YLM{^)TLiz;q,p)YlMip), (3.33) 

LM 

where q = |q|, p = |p| and q and p denote unit vectors in the direction of q and p, 
respectively, and YiMi^l) = (47r)^^/^ (q|LM). Since the ground state contains only S-wave, 
only the term with L = counts in these sums. Hereafter, we shall suppress the index L in 
the partial wave amplitudes. 
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Below we assume, that potential ^(q, p) is short-ranged, hermitian and real. In this case, 
the unitarity condition for the scattering matrix t{z) in the S'-wave then takes the form 

lmt{z;q,p) = -^t{z;q,ko)t*{z;ko,p), ko = ^ME . (3.34) 

Now, solving Eq. ( ^.321 ) by iterations, for the energy shift we obtain 

= Re (^olt(^o)l^o) + Re (*o|t(^o) g'"''(^o) t(^o)l^o) 

+ Re {^o\t{Eo) g^'-^iEo) t{Eo)\^o) + «(«') ■ (3.35) 

Below, we proceed with the evaluation of all three terms in the right-hand side of 
Eq. ( p35D : 



1. The term with no Coulomb exchanges is given by 
Re(^o|t(^o)|^o) = 0oRei(^o;O,O) 

+200 / rfz/(q)^o(q) (Re t(^o; 0) - Re t{Eo; 0, 0)) + o(a^) , (3.36) 

where 0o = / diyi^)'^oi^) = il^/'^Y^'^- Using the following property of the scattering matrix, 
that can be proven for short-range potentials (see Appendix ^ 

lim- (Ret(Eo;q,0) -Ret(EQ;0,0)) = 0, (3.37) 

q^O q 

this term can be rewritten as 

Re(^o|t(^o)|^o) = 00 {Ret(^o;0,0) + ^ Q[t]} + o(a^) . (3.38) 
Here t = Ret(2M; 0, 0) and 

Q[t]=J^{Ret{2M;q,0)-t) . (3.39) 







Note that we have substituted Eq = 2M in the correction term, since this does not affect 
the result at the accuracy we are working. 

Next, we have to perform an analytic continuation of the scattering matrix from Eq to 
2M. The distance between these two points is of order a^. However, since the real part of 
the scattering matrix has the unitary cusp at threshold, the difference between t(-Eo) and 
t{2M) is of order a rather than a^. It is easy to check that 

aM^ - 

Re ^(£"0; 0, 0) - t = t 2 + o(a) (3.40) 

87r 

Collecting all contributions, we finally get 

Re(^o|t(i?o)|^o) = <Pl\t + -^t 2 + Q[t] + o{a^) . (3.41) 
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2. In the calculation of the matrix element corresponding to 1 Coulomb photon exchange, 
it is convenient to separate the contributions coming from small and large integration mo- 
menta. This can be achieved by rewriting the Coulomb potential as 



1 



1 



1 



1 



|q — pP |q — pP + 6^ J 



|q — pp |q — pP + 6^ 
The integration can be straightforwardly carried out, resulting in 



(3.42) 



(^ro|t(Eo)g'-''(^o)t(£^o)|^o)=0o<j 
where 

^ dpdk {p — k)'^ + l? 



aM^-^^ b 



4:71 



+ 



27 Att^ 



n[b; t] \ + o{a' 



(3.43) 



n[b; t] 



+ 



In- 

pk [p + kY + b"^ 
dpdk 



t(2M;0,p)t(2M; k, 0) 



(3.44) 



pk 



1 (P - kY , ip - kf + 6^ 
In -, —- — m ■ 



{p + kf 



{p + ky + 62 



t{2M;0,p)t{2M; k, 0) - t 



It can be checked that the result does not depend on the arbitrary cutoff parameter b. 
3. The integration in the matrix element containing many-Coulomb Green function, can 
be directly carried out. The result is 



3aM^ 



6TT 



Finally, putting things together, for the energy shift we obtain 



AE 



An 



Air 



t Mn h 

27 



AaM , aM^ J , ,^ 



(3.45) 



(3.46) 



The functionals 7^ [6; t\ and Q\t\ are given above. The expression for the energy shift can be 
rewritten in form similar to Eq. ( p.25|) 



AE 



a' 



A 



NR 



in 



An 



^^^(ln«-l)^ + o(«^ 



(3.47) 



where 



A 



NR 



aM 



8. — SW-l^'^lMl+oH. 



(3.48) 



In order to check the universality, one has to evaluate the full scattering matrix T which 
is defined by Eq. (|3.29|) , in the vicinity of threshold. We use the scattering theory on two 
potentials, and evaluate T perturbatively at 0{a) 



T = hc + t + hcgot + tgohc + tgohcgot + o{a) , 



T = T-h 



c ■ 



(3.49) 
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The integrals entering here, are calculated similar to given above. Note that, even 
there are no more ultraviolet divergences, we still use dimensional regularization in order 
to regularize infrared divergences. After removal the Coulomb phase, one may safely put 
(i — > 3. The truncated amplitude in the vicinity of threshold behaves as 



Re[e 



-2iaec 



;q|T|p)] 



£;=q2/A/, |q| = |p|->0 



4:71 



Pin 



2M 

M 



^^^ + 0(|p|), (3.50) 



with exactly the same A^^ as in Eq. (|3.48|) . This means that, we have verified the univer- 
sality conjecture formulated in the beginning of this section. 



C. Matching condition for the potential 

The universality conjecture, proven in the previous section for the one-channel case, 
provides one with the matching condition for the short-range potential M(q, p). In order to 
derive this condition, we assume that short-range potential, in analogy with Eq. ( 3.26| ), can 
be written as 

u{q, p) = Mo(q, p) + aui{q, p) + O(a^) . (3.51) 

Further, we introduce 

to = Uo + Uogoto. (3.52) 

Using scattering theory on two potentials, we obtain 

t = to + tt (1 + ^ogo) ui (1 + goto) + 0{a'^) = to + ati + O(a^) . (3.53) 

The matching condition for the potential can be obtained by using Eqs. p.23| ), ( |3.26D and 
( p.48| ) order by order in a 

A^ = -4MHo, (3.54) 

= -4M2 |(1 + togo)^/i(l + goto) - ^ t^ In — + — Q[to] + ^ 7^[6; to] j , 

where to = Reto(2M; 0; 0). 

It is seen, that the matching condition imposes rather loose constraints on the potential 
M(q, p). For example, the matching condition in the first line requires that the strong 
scattering lengths in the relativistic theory and in the potential model are the same. The 
behavior of the scattering matrix above threshold does not play any role. This property of the 
matching condition is not surprising, if one adopts the interpretation of the potential model 
given in the previous section, namely that the short-range potentials are regularizations of 
the contact interactions that arise from field theory. The looseness of the matching condition 
then merely reflects the freedom in the choice of such a regularization. Using this freedom, 
one may further specify the potential, assuming 

«i(q,p) = Auo(q,p) . (3.55) 
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This amounts to the matching of the "strength" of the potential, that is the counterpart of 
the couphng in the non-relativistic Lagrangian, at 0{a), with the momentum dependence of 
the potential fixed by hand. Using the above ansatz, one can rewrite the matching condition 
in the following form 

A? = -4M^ {a (f„ - ^ew) - + ^ 21'°) + j^^'^l''-*"!} ■ (3.56) 

where 

g[t]= dpt{2M;0,p)t{2M;p,0) . (3.57) 
Jo 

This matching condition uniquely determines the couplings in the potential at 0(1) and 
0{a). Namely, first, one chooses the momentum dependence of the potential Mo(q, p) and 
adjusts the coupling at 0(1), to reproduce the strong scattering length obtained from the 
relativistic theory (the first line of Eq. ( p.56|) ). This completely determines the scattering 
matrix at all momenta and energies. At the next step, with a given t^, one calculates all 
integrals entering the matching condition at 0(a) (second line of Eq. ( p.56|) ), and determines 



the coupling constant A from this equation. Despite the unrestricted freedom in the choice 
of the momentum dependence of 'Uo(q, p), the potential M(q, p) that we construct, repro- 
duces, by construction, the strong energy-level shift at 0{a^), as well as threshold scattering 
amplitude at 0(a). 



D. Separable potential 

For demonstrative purposes, below we again consider the solution of matching equation 
for simple rank-1 separable potential 

Uo = gu{q)u{p)] u{q) = -— — - (3.58) 

With this potential, one is able to carry out the calculation of all quantities entering the 
matching condition, in a closed form. The result is given by 

n[b; to] = - 17 In , g[to] = ^Pto' (3.59) 

Substituting these expressions into the Eq. (|3.56|) , one obtains two equations which fix 
the couplings g and A at 0(1) and 0(a), respectively 



4M2 1 + 



32ttM 



1 [A? 2M 1 p , , ^ 



327rM 
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As expected, the potential range parameter (3 is not determined from matching condition. 

We mention again, that the separable form of the strong potential ( 3.58 ) is not, by far, 
the unique choice. E.g., the local square- well strong potentials used in Ref. [0], can serve 
equally well. In the latter case, however, we were not able to obtain the corresponding 
integrals that enter the matching condition, in a closed form. This means, that numerical 
methods should be used. 



IV. TWO-CHANNEL CASE: DECAY OF THE PIONIUM 

In the previous Section we have considered the strong energy shift of the hadronic bound 
state at 0{a^). The problem was solved in the one-channel model - that is, the theory 
contained, from the beginning, only a doublet of charged scalar particles. In addition, we 
have neglected relativistic corrections and higher-order derivative couplings in the strong 
Lagrangian. As was mentioned, all these effects do not contribute to the strong energy shift 
at the accuracy considered, and the (possible) coupling to other channels does not show 
up explicitly in the energy-shift calculations at 0{a^) as well pTI]. The situation is quite 



different, if one considers the decay of the hadronic bound state at next-to- leading order in 
isospin breaking: here all above effects should be consistently taken into account. In this 
Section, we shall enlarge the formalism developed in the Section |T| to the two-channel case, 
and include relativistic effects as well as derivative interactions. To this end, we consider 
the decay of the 7r~^7r~ atom (pionium) which has been already investigated in detail within 
the non- relativistic effective Lagrangian approach and ChPT |ll9|j20| , |22[| . Generalization to 



other hadronic atoms, as well as inclusion of spin effects is straightforward an will not be 
discussed. 

Preliminary remarks are in order. The isospin breaking in the tt^tt^ system is due 
to two physically distinct sources: electromagnetic corrections which are parameterized by 
the fine structure constant a, and the quark mass difference rrid — rriu. The notions of 
"isospin-symmetric world" and "pure strong interactions" refer to the idealized world with 
a = 0, — rriu = 0, and where, by convention, the mass of the pion coincides with the 
charged pion mass in the real world. It is convenient to introduce the common counting 
for two different isospin breaking parameters. In the following, we shall use the following 
assignment p9|j20| j2^ 



a^{md-mu) (4.1) 

The TT+TT^ atom decays predominantly into the 7r°7r° final state: Ttot = ^2itO + 
0{S^) ||T9|j20|j2^ . At leading order, r27rO = 0(5''/^). We shall be interested in next-to- leading 
order corrections in isospin breaking - up to and including order 5^/^. At this accuracy, 
still only the decay into Tr^vr" final state is possible [Q, so it is perfectly consistent to re- 
strict ourselves to the consideration of the two-channel (tt+tt" and tt'^tt'') problem in the 
quantum-mechanical framework. 

Our strategy will be follows. First, we consider the "relativization" of the non- relativistic 
effective Lagrangian approach used in Refs. |T^,pD|,p2[ , in order to bring this approach in con- 



formity with the relativized potential model used in Refs. 0-0 to treat the same problem. 
The relativized field-theoretical framework is used to calculate the decay width of the 7r"'"7r~ 
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atom at order 5^/^ - the result, of course, agrees with that from Ref. |jT9|- Next, we con- 
struct the two-channel potential model which reproduces the result of the field-theoretical 
approach. The universality conjecture is verified for the two-channel potential model, with 
relativistic corrections and derivative interactions taken into account. This provides us with 
the matching conditions (one per channel) for the short-range potential, which can be solved 
similar to the one-channel case. 



A. Relativistic corrections 



In Ref. [£2| it has been argued that the following non-relativistic effective Lagrangian is 



sufficient to carry out the calculation of the decay width of the tt+tt atom at 0(5^/^) 
C = Cq + Lb + Cc + Cs 

A 



^z.= E + (4.2) 



Cc = -47r«(7rl7r_)A"^(7r;'^7r+) H , 

Cs = cinlnln+n^ + C2[7r^7rl(7ro)^ + h.c] + C3 (tTj^ttq)^ 

tt^ttItto a ttq + h.c] + ■ ■ ■ . 

This Lagrangian is built from the non-relativistic pion field it ±,7^0- It contains the non- 
relativistic kinetic term Co, along with the term Cd which accounts for the relativistic 
corrections to the pion energy. These corrections have been included in the bound-state 



width in a perturbative manner [1T9| , P2| . Further, the formation of the bound state proceeds 
mainly due to the static Coulomb interaction contained in Cc, whereas strong interactions 
described by the local four-pion Lagrangian Cs are mainly responsible for its decay. The 
constants in Cs are determined from matching to the relativistic theory. Again, as in the 
one-channel case, we have truncated all terms which do not contribute to the quantity of 
interest (the decay width) at the accuracy we are working. 

In order to bring our framework in conformity with the relativized potential model which 
have been used for the study of pionium decay [^|-|^ , below we include the relativistic correc- 
tions contained in Cd, in the unperturbed Lagrangian. Diagrammatically, this corresponds 
to summing up all mass insertions in the free non-relativistic pion propagator (see Fig. In 
actual calculations of the diagrams in effective field theory, one is again forced to treat part 
of these corrections perturbatively, in the expansion of powers of p^/M^,. Since our trick 
amounts merely to the redistribution of terms in the total Lagrangian between unperturbed 
and interaction pieces, the results for any observable quantity (e.g. the decay width) should 
remain unaffected. The reason why this redistribution is carried out, is twofold. 

i) If one has the same bound-state equations in the effective field theory and in the 
potential model, one can merely read off the potential from field-theoretical bound- 
state equations. 
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ii) Though the perturbative treatment of the mass insertions in the effective field theory is 
easy, this becomes rather comphcated in the potential model with general non-contact 
interactions. Technically, it is preferable to have a framework where these insertions 
are summed up from the beginning. 

In order to design such a framework, we bring together Cq and £d to form the relativized 
kinetic term 



Ci, = Co + Cd= E 4 Udt - vr, . (4.3) 

i=±,0 ^ ^ 

The corresponding free relativized Hamiltonian is given by 

= / ci^x 4(0' V^'- A 7r,(0,x) , (4.4) 

and Coulomb (He) and strong (H5) Hamiltonians are defined in analogy with Eq. (|3.4| ). 
The full scattering matrix obeys the Lippmann-Schwinger equation 

T\z) = (Uc + Us) + (He + Us) G'iz) T'{z) , G'(^) = . (4.5) 

Z — -Hr 

Poles of the scattering matrix on the second Riemann sheet of the complex z-plane 
correspond to unstable bound states. The real and imaginary parts of the pole position 
determine the energy and width of such a bound state, according to E = Iiez,T = —2lmz. 

Further, we define the transformed quantities 

, XT \ 1/2 



.(Q,q|T(--)|P.p). ^ .(Q,q| (i±^) T'(.) {i±^) |P,p), 

.(Q,q|V-W|P,p).= ,(Q,q|(i±^)"^Hc + H.)(£±|^)"^P,p), (4.6) 



^(Q,q|Go(^)|P,p)B= A(Q,q| . /I A G\z) . J |P,p)b. 

\sa[z + \iiR) ) \sb[z + \±r) ) 

with s+ = 1, So = 2! and A, 5 = +, label the charged (tt+tt^) and neutral (vr^vr'^) channels. 
In the CM frame, the free resolvent (see Eq. ( [4.6|) ) takes the form 



A(q|go(^)|p)B = (27r)353(q_p 



\ 

.f^' , 

V - / AB 



(4.7) 



where 



We see that, despite the non-relativistic appearance of the free resolvent ( [4.7|) , the relativistic 
effects are taken into account in the definition of the quantities fc^, fcg, Eq. ( [4.8|) . 



kl='^-Ml,^ kl='^-Ml,. (4.8) 
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The Lippmann-Schwinger equation for the transformed scattering matrix now exactly 
corresponds to the Schrodinger equation used in the relativized potential model 



i{z) = v{z) + v{z)go{z)i{z) . (4.9) 

If one expands the transformed potential in Eq. ( [4 .61) in the small quantities p^/M^^, 
z — 2M^+ ~ 0(5^) and z — 2Mt^o ~ 0(5), it is seen that the transformation formally 
amounts to the following replacements 

r z - 2M^+ z - 2M^o 1 

c.^c.-^ + .... (4.10) 

Since the matching condition determines these couplings from the amplitudes evaluated 
at threshold, one has to substitute z = 2Mj,+ or z = 2M^o (depending on the particular chan- 
nel) in the above expressions. Finally, one may conclude that the transformation, introduced 
above, amounts to the following redefinition of strong couplings in the Lagrangian ([4.2|) 



(4.11) 

Accordingly, one may define the energy-independent potential V which is obtained from 
the transformed potential V{z) by substituting the threshold values for the parameter z. In 
the CM frame this potential is given by 



Ana 



I -c'l v^(-4 + 2cUq2 + p2))' 

(q|^|p) = - i "^'"12 + (4-12) 

1^ PI \^V2(-4 + 2cl(q2 + p2)) -24 



By construction, the energy-independent potential when used in the Lippmann-Schwinger 
equation ( [4.9| ), reproduces the decay width in the first non-leading order in isospin breaking. 

Note that the relativization scheme described above is linked to the particular choice 
of the Lippmann-Schwinger equation in Ref. If a different choice is assumed, it is 
straightforward to adapt the above relativization scheme to this new equation. 



B. Decay width of the pionium 

One may study the pionium decay in the relativized framework by using exactly the 
same methods as within the non-relativistic effective Lagrangian approach. We refer reader 
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to Refs. |lT9| , pO| , |22[] for more details concerning the technique, and merely quote the final 
results here. The decay width is given by 

T2^o = laYA'{l + K). (4.13) 

where p* = {M^+ — M^o — | a^M^+)^/^ and the quantities A and K are defined below. 
The matching condition now reads 

Aici\iiz)\p)B = - / ^ ^ ^ TU^, p) ^ , A,B = +,0. (4.14) 

2^MASA(^A{q) 2^MbSbUJb{P) 

The threshold behavior of the relativistic tt^tt^ n^n^ scattering amplitude at 0{6) is 
given by (cf with Eq. ( p.22 )) 



Re(e-^°''^T^o(q,P)) 



' 5f ln^ + Re<7°° + 0(|q|). (4.15) 



IqHo |q| 

And the quantity A is, from the matching condition, expressed through the regular part of 
the relativistic tt+tt" scattering amplitude at threshold 

A = Re A+7°° + o{6) , A = ao-a2 + 0{6) , (4.16) 

where a/ denote the S'-wave vrvr scattering lengths in the isospin limit in the channel with total 
isospin /. Further, the quantity K is expressed through these scattering lengths according 
to 

K=-{ao + 2a2f - — (Ina - l)(2ao + a^) + o{6) . (4.17) 

with K = Ml+/Mlo - 1. 

The matching condition for the couplings c'^^ ■ ■ ■ c'^ is given by 



An 
3M! 



t:7I 

c'l = ^772- (2«o + 02) , 



4 = («o + 2a2) . (4.18) 

This matching condition determines c[,C2,c'^ at 0(1). In addition, it fixes a particular 
linear combination of c'2 at 0{6) and C4 at 0(1). Having expressed these couplings through 
the relativistic scattering amplitude, one may again interpret the short-range part of the 
energy-independent potential (|4.12|) as the contact potential to be used in the Schrodinger 
equation - provided one uses the dimensional regularization throughout. However, in a 
complete analogy with the one-channel case, there is also a possibility to construct the 
conventional short-range potential which reproduces the field-theoretical result for the width 
at 0(5^/^) and does not lead to the ultraviolet divergence. This is demonstrated in the 
following Section. 
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C. Pionium decay within the potential model 



In this Section, we shall demonstrate that universality, which we have proven for the 
ground-state energy-level shift in the one-channel case, also holds for the decay width of the 
pionium. Namely, we prove that the relation between the decay width of the ground state 
and the threshold amplitudes is the same in field theory and in the potential model in the 
first non-leading order in isospin breaking. Given the fact, that the energy and the decay 
width are the only observable characteristics of the bound state, one may conclude that the 
properties of the bound state in the first non-leading order in isospin breaking are completely 
determined by the properties of the amplitudes. 

The calculations which are performed below, are analogous to the ones already carried 
out in the one-channel case (see Section |T|). We give only the final results here. 

The scattering matrix in the CM frame obeys the Lippmann-Schwinger equation 

i{z)=v + vgoiz)i{z). (4.19) 

All quantities entering this equation are 2x2 matrices with the entries A, B = +, 0. In 
particular, the potential of the above Lippmann-Schwinger equation is given by 

Ai<i\v\p)B = 5^+5^+ + «^^(q, P) = (q|hc|p) + u^^(q, p) , (4.20) 

|q- p| 

where ^^^(q, p) stands for the short-range part of the potential. 

The two-channel counterpart of the scattering matrix t{z) introduced in Eq. ( 3.13 ) sat- 
isfies the following equation 

t{z) =u + ugc{z)T{z) , (4.21) 
where the pole removed Coulomb Green function is given by 

nq\gc{'2M^+ + kl{z)/M^+)\p) 

A{q\gc{z)\p)B = 

V (q|go(2M,o + A;2(z)/M^o)|p) 

Further, in analogy with Eq. ( p.31| ), we again introduce 

t{z) =u + u go{z) t{z) , t{z) = t{z) + t{z) 6g{z) t{z) , 6g{z) = gc{z) - go{z) . (4.23) 

We again, as in the one-channel case, assume, that the potential is hermitian and 
real. In addition, we assume that u"^^ = u^^. The unitarity condition for the scattering 
matrix t{z) in the two-channel case takes the form 

Imt^^(z; g,p) = -^^^^^ t^^z; q, k^{z)) {t^^{z- k^{z),p)y 

_ Mi:pi^t^o^,,q^k,{z)){t'^{z-k,{z),p)y, (4.24) 
An 

and the partial- wave expansion is again performed, according to formulae (|3.33| ). 



(4.22) 
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1. Decay width 

The position of the bound-state pole in the full scattering matrix is determined by the 
equation (cf with Eq. ( p.l2| )) 

z-Eo-(^o|r++(z)|^o) = 0. (4.25) 

Solving this equation by iterations, it is straightforward to demonstrate (see Appendix jB]) 
that the decay width F = —2\m.z can be rewritten in a form ([4.13|) , if we identify 



^ = — 4v^ m":mH' I Re i +° - ^ In Re t +°Re i 
327r ^ \ An 

+ Q[t] + t] I + 0(5^) , (4.26) 

where the following notations are introduced 

t++ = t++(^+; 0, 0) , r° = t+°(;.+ ; 0, A:o(^+)) , P = t°°(^o; 0, 0) , (4.27) 
and z+ = 2M^+ and Zq = 2M^o. The functionals Q[t] and TZlb; t] are defined in Appendix 0. 

5. Matching condition for the threshold scattering amplitudes 

In order to check the universality conjecture, we calculate the scattering amplitude for the 
process tt~^tt~ vr^vr^ in the vicinity of threshold, in the potential model. The perturbative 
expansion in a for the scattering matrix has the form 

i +° = t+° + he (go)++ + t++ (go)++ he (go)^"- + o{5) . (4.28) 

In order to calculate the integrals entering here, one has to use the dimensional regulariza- 
tion to regulate the infrared divergence caused by the one Coulomb photon exchange. At 
threshold, the behavior of the scattering matrix is given by 



Re[e-^"^^ +(q|t|p)o] 



""^'^^t+^-^Ret^+Ret-'-o In-''^ 



|qHo,|pHfco(^+) 4|q| 4:71 

1 327r 

+ r- 3/2 1/2 ^ ^ + 0(|q|) . (4.29) 
4:^/2 mI'+ My 3 



In addition, we shall use the matching condition in the isospin limit in the elastic channels 

TT + TT^ TT + TT^ and 71^71^ Tl^Tl^ 

4:71 - 4tt 

Ret ++ = (a2 + 2ao) + 0(5) , Ret °° = (ao + 2a2) + 0{5) . (4.30) 



One may straightforwardly check, that the regular term A in Eq. ( 4.29|) is indeed given by 



the expression ( |4.26| ) obtained from bound-state calculations. Consequently, the universality 



conjecture holds also in the two-channel model, in the context of the bound-state decay 
problem. 
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D. Matching procedure for the potential 



Using the exphcit expression (|4.26|) , we shall derive the matching condition for the short- 
range potential by matching the amplitudes at threshold in the relativistic theory and in 
the potential model. The universality then ensures, that the bound state decay width is the 
same in both theories at 0(5^/^). 

As in the one-channel case, we assume that the short-range potential contains the isospin- 
conserving and isospin-breaking parts 

«^^(q, p) = «^^(q, p) + <^(q, p) + 0(6') , ~ 0(1) , <^ ~ 0(6) , (4.31) 

and the matching condition holds separately at 0(1) and 0{6). 
It is convenient to introduce the purely strong scattering matrix 

to{z) =uo + uogo{z)to{z) , (4.32) 

where go denotes the free resolvent in the absence of isospin breaking 

A(q|go(^)|p)B = (2vr)=^53(q - p) 5^^ f^^^ ^ . (4.33) 

The perturbative solution for the scattering matrix t in terms of to is given by 

t = to + (1 + togo) ui (1 + goto) + to(go - go)to + to(go - go)to(go - go)to H , (4.34) 

where ellipses stand for the terms that do not contribute to the real part of the scattering 
matrix t at 0{5). One has to mention here, that the corrections in Eq. (|4.34| ) containing 
go ~ go 5 are usually referred to as "mass-splitting corrections" within the potential model. 
These corrections are, of course, absent in the one-channel case. Moreover, it is implied in 
the existing potential models, that these corrections fully account for the effects caused by 
the charged and neutral pion mass difference, that amounts to the negligence of the mass 
corrections contained in the relativistic amplitude. 

On the other hand, the relativistic scattering amplitude that enters the matching condi- 
tion, can be written as p0| , |22[| 



^ = ag - 02 + e + 0(52) ^ (^4 35^ 

where, in order to calculate e one invokes, e.g., ChPT. 

Substituting these expressions into the matching condition for the amplitude in the 
charge-exchange channel, we obtain 



327r 1 



Hh] , (4.36) 



where ^[to] is the complicated functional depending on to- Its explicit form is given in 
Appendix y. 

The matching condition in the charge-exchange channel should be complemented by two 
additional conditions from elastic channels, where it suffices to perform matching at 0(1). 
Finally, at 0(1) the matching condition reads 
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4y2vr \ _-+o 

TV' 

(2ao + aa) = C+ 



.T + 

An 



(ao + 2a2) = to (4-37) 



3M2 

At order 5, from Eq. (|4.36| ) we obtain the following matching condition 



'e + -(ao-a2) = {(1 + to go) t^i (1 + go ^o)}"^ 



3M\ V 4 



+ ^[to] . (4.38) 

thr 



As expected, the matching condition does not pose any constraint on the momentum 
dependence of the short-range potential. We may use this freedom and assume 



<^(q,p) = ^<^(q,p). (4.39) 

With this assumption, the integral entering the matching condition at 0(5), can be 
expressed through to 



thr 



-UV°-^^[^o]L (4.40) 



27r2 



{(l + togo)«i(l + goto)}+° 
where 

oo 

^[to] = jdp {t++(^+;0,p)t+°(-2+;i',0) + t+°(2:+;0,ij)C(^+;P,0)} , (4.41) 



and the matching condition then uniquely determines the coupling uj in the short-range 
potential. 

The matching procedure proceeds in several steps 



i) One solves the equation ( 4.32|) in the absence of isospin breaking, in the basis of isospin 



scattering matrices to, t2, and adjusts the couplings in the short-range potentials so, 
that the matching conditions ( [4.37| ) at 0(1) are satisfied. The scattering matrices in 
physical channels are expressed via to, t2 according to 

= I (2to + 12) , C = ^ - to) , t°° = 1 (to + 2t2) . (4.42) 

ii) With the given t^^, t^^, t^^ one calculates the functional JF[to] that appears in the 
right-hand side of the matching condition at 0{6) ( [4.38| ). If the potential in each 



channel is given, e.g. in the simple separable parameterization ( 3.58 ), this again can 
be done in a closed form. For the general short-range potential, one, however, may 
have to resort to the numerical methods. 

iii) Finally, one determines the coupling uj in the short-range potential from the solution 
of the matching equation at 0{6). 
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The two-channel hadronic potential that is obtained in a result of the matching condition 

reproduces, by construction, the relativistic amplitudes at threshold. The 7r"'"'7r~ 7r°7r° 

amplitude is reproduced at 0{S), whereas vr+vr" tt+tt", tt^tt^ — 7r°7r° at 0(1) (no isospin 

breaking). Moreover, the potential reproduces the decay width of the tt+tt^ ground state at 
0(59/2)_ 

Note that, if one is willing to reproduce amplitudes in all three channels at 0{5) at 
threshold, as well as to describe simultaneously the energy shift and the decay of the tt+tt" 
bound state, the simple ansatz ( [4.39|) does not suffice. The possible generalization of this 
ansatz is given by 

Ui = Ouo + Uod, O = CJil + 6^20"! + U^sfTs , (4.43) 

where 1 denotes the unit 2x2 matrix, cTj are Pauli matrices, and real parameters ui, 0^2, uj^ 
are the couplings in the short-range potential at 0{S). The generalization of the matching 
condition at 0{6) to the elastic channels with the use of the above ansatz is straightforward. 



V. COMPARISON TO THE EXISTING POTENTIAL MODELS 

As was discussed in the Introduction, the calculations of the observables of hadronic 
atoms carried out within the potential approach (see e.g. ||7HlO|) often lead to the results 



which are in a pronounced disagreement with those obtained on the basis of ChPT. In 
addition, these models are sensitive to the particular choice of the interaction potential, 
and/or to the choice of the unperturbed Green function in the Lippmann-Schwinger equation 
- the different choice may sometimes lead to the dramatic consequences . The reason for 
this is now clear: the existing calculations use the potential which is not matched to ChPT 
in the isospin-breaking phase. Instead, e.g. in Refs. the strong potential is matched 
to ChPT phase shifts up to 500 MeV that, in our terminology, corresponds to the matching 
of the parameters of the effective-range expansion in the strong sector at a high order (see 
Section |lD. In the contrary, the isospin-breaking part of the short-range potential is put 
to zero by hand: it is assumed that the entire isospin-breaking effect in QCD is due to the 
Coulomb interactions, kinematical effects due to the mass difference between charged and 
neutral particles, etc. However, it is clear from our construction that such an assumption 
already at the next-to-leading order ignores some contributions which are present in QCD 
and are contained in the isospin-breaking part of the amplitude (Eq. ( |3.26| )). Below, we 
list several of these contributions. 

1. Scattering amplitudes in QCD explicitly depend on quark masses - this dependence 
is governed by the chiral symmetry. For demonstrative reason, let us consider the tree-level 



scattering amplitude ir'^n — > ir^iT^ in ChPT (see Ref. [22| for more details) 



rjiR -5 — 2rhB 

where s is the Mandelstam variable, m = |(mu + m^), and the constants F and B are 
related to the pion decay constant and the quark condensate in the chiral limit. At order 
in ChPT, 2mB = M^o- According to our convention, however, the pion mass in the isospin- 
symmetric world is equal to M^+, and the splitting of the above scattering amplitude at 
threshold into the isospin-conserving and isospin-breaking parts takes the form 



28 



where the isospin-breaking correction (second term of Eq. ( |5.2| )) arises due to the exphcit 
dependence of the scattering amphtude (|5.1|) on rh. Furthermore, the relativistic Lagrangian 
that leads to the scattering amphtude ( |5.1| ) in the tree approximation, should also contain 
m 

= ——p^ if'^ ip~ (p^ ip^ + terms with derivatives . (5.3) 

On the language of the potential theory, this would correspond to the short-range potential 
with the coupling proportional to rh. This coupling changes when one goes from the isospin- 
symmetric world with 2m' B = to the real world 2mi? = M^o (the constant B is chosen 
to be the same). Put differently, in order to take the above effect into account, the short- 
range potential should contain the isospin-breaking piece. Since it is not the case in the 
existing potential models, this effect is missing there. 

2. Direct quark-photon interactions in QCD that occur at a QCD scale (around 1 GeV), 
correspond to local vertices in ChPT (see Fig. H). The typical four-pion Lagrangian has the 
form similar to (|5.3| ) 

= e^ki {p~^(p~ip^ip^ + terms with derivatives , (5.4) 

where ki are the so-called electromagnetic LECs. The potential which corresponds to this 
interaction, obviously vanishes in the isospin limit e = - consequently, it is not included in 
the existing potential models as well. 



3. Virtual photon corrections like ones indicated in Ref. ||25| may lead to the serious 
discrepancy with the predictions of the potential model. The counterpart of such a diagram 
for TTTT scattering process is depicted in Fig. ^ It is clear that, at low energy, the contribution 
of this diagram is regular in kinematical variables. Consequently, it should be included in the 
isospin-breaking part of the short-range potential. Since the existing potentials are assumed 
to be isospin-symmetric, they miss this particular contribution as well. 

There is a number of other effects which are not properly included in the potential 
models (e.g. the kinematical mass-splitting effect in t and u channels), as well as the effects 
which are correctly treated (e.g. resummation of Coulomb ladders). As a general rule, the 
potential model gives a reliable prediction of the isospin-breaking effects if and only if the 
potential is matched to ChPT in the isospin-breaking phase. We note that the effects that 
were discussed in this Section, are not by all means small. In fact, they form the bulk of the 
total isospin-breaking corrections both in the tt+tt" atom width and the 7r~p atom energy 
shift. It is worth to mention that, in addition to the isospin-breaking effects, the matching 
takes care of the differences caused by choice of the potential and of the free Green function. 
In particular, there is no necessity to fit the strong phases up to high energies: it suffices to 
fit only the scattering lengths. 

From the discussion above, it is clear that, strictly within its range of applicability, the 
potential approach does not provide us with a new physical information as compared to 
ChPT. In the view of the extremely complicated analysis of the isospin-breaking corrections 
in the scattering processes [|1^], it can be still useful, for the time being, to construct the po- 
tential with full isospin-breaking content of QCD at threshold. In the absence of systematic 
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analysis carried out on the basis of ChPT, one may hope that this construction would allow 
one to improve the quality of predictions of isospin breaking for the scattering experiments. 



VI. SUMMARY 

i) On several examples with an increasing level of complexity, we have demonstrated the 
guidelines for the derivation of the short-range hadronic potential from underlying field 
theory of strong and electromagnetic interactions. These potentials, when inserted 
in the conventional Lippmann-Schwinger equation, reproduce, by construction, the 
threshold scattering amplitudes and observables of hadronic bound states in the first 
non- leading order in isospin breaking parameter (s). 

ii) According to the viewpoint adopted in the present paper, the conventional short-range 
potentials are considered as a mere regularization of the singular pointlike interactions 
which describe low-energy interactions of hadrons in the field theory. For this reason, 
the shape of the potential does not bear a physical information. Couplings in the 
potential are determined from matching of the scattering amplitudes in the potential 
approach and in the underlying field theory, both expended in powers of the CM 
momentum squared p^. Performing the matching at higher order in p^, one arrives 
at a more accurate description of both the scattering amplitudes and bound-state 
observables in the potential approach. Apart from the above matching condition, no 
further restriction is imposed on the potentials. 

iii) The derivation of the potential for the description of bound states is based on the 
universality, which states that the properties of bound states are the same in the 
potential scattering approach and field theory, once the threshold scattering amplitudes 
are the same. Due to the universality, one may carry out the matching in the scattering 
sector, where the perturbation expansion in a works. 

iv) The reason why the results of calculations for the observables of hadronic atoms carried 
out within the potential approach ||3"10| generally differ from those obtained in ChPT, 



is now crystal clear. In order to agree with the latter, the potentials should be matched 
to QCD in the isospin-breaking phase. The matching, in general, generates a nonzero 
isospin-breaking part of the short-range hadronic potential. Furthermore, it may turn 
out that the prediction for the isospin-breaking corrections to the atom observables 
in the potential approach is close to that of ChPT 0. This means that the isospin- 
breaking part of the short-range potential (only in this particular hadronic channel) 
constructed through the matching procedure, is very small. 

It is obvious that, in the context of hadronic atoms, the potential which we have 
constructed does not contain new physical information as compared to already avail- 
able solution of the problem on the basis of effective chiral Lagrangians. However, 
it is interesting to study, whether the constraints imposed on the potential from the 
matching of the amplitudes at threshold, can be useful in the analysis of the scattering 



data along the lines similar to those from Refs. 
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APPENDIX A: BEHAVIOR OF THE SCATTERING MATRIX AT A SMALL 

MOMENTA 

In this Appendix, we consider the behavior of the scattering matrix at a small momenta 
in the potential scattering theory with a short-range potential u. The S'-wave potential in 
momentum space is given via the Fourier transform 



ridri r2dr2 u(ri,r2) , (Al) 

J a p 



where u{ri,r2) is assumed to be short-ranged. The expansion of the right-hand side of 
Eq. (|A1|) in powers of q contains only even powers of q, provided the integrals in ri which 
emerge after the expansion of the integrand, are convergent. In particular, 

lim ^(^^P)-^(^^P) _ Svr^ 7^4^^ 7^ sin(pr2) 



/ dri / r2 dr2 u[ri, r2) , (A2) 

^ 



q^^o g2 3 J ' J p 

if the integral in the right-hand side is convergent. We assume that our short-range potential 
satisfies this requirement. 

The scattering matrix satisfies the Lippmann-Schwinger equation 

t{E; g, p) = u{q,p) + —l ^ _ ^^^^ u{q, k) t{E; k, p) . (A3) 

Using this equation, it is immediately seen that the scattering matrix has the similar 
behavior at small g, as the potential 



,^,m,,P)-mM (A4) 



The equation ( |3.37| ) directly follows from this equation. 
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APPENDIX B: DECAY WIDTH IN THE TWO-CHANNEL SCATTERING CASE 



The perturbative solution of Eq. (|4.25|) gives 

z-E^- (*o| t++(£;o) l^o) - (^ol t++(^o) g'-^'l^o) t++(^o) l^o) 

- (^ol t++(^o) g^^^'l^o) t++(^o) l^o) + ■ ■ ■ = . (Bl) 

The separate contributions in this equation correspond to sphtting of the Coulomb Green 
function into the no-photon, one-photon, and many-photon pieces. 

In the evaluation of the matrix element with no photons, we are again confronted with 
the necessity of the analytic continuation of the scattering matrix from the bound-state 
energy Eq to the scattering threshold z+. Due to the presence of the unitary cusp in the 
real part of the scattering matrix, in a result of such a continuation we get the difference 
which is non- vanishing at 0{5) 

Ret+°(Eo; 0, A;o(Eo)) - Ret+° = t+o + o{5) . (B2) 

OTT 

The imaginary part of the matrix element with no Coulomb photon exchange is given by 
Im (vl/ol t--{EQ) \^q) = ||,-.0|. 

47r TT 



where 



Q[t] = j% (Ret+°(;.+ ;p, ko{z+)) - Ret+°) . (B4) 
^ 

The imaginary part of the matrix element with one Coulomb photon exchange is given 



by 

Im(^o| t++(^o) g'^^'l^o) t++(^o) 1^0) = 
47r I 271 27 J 



(B5) 



where 



X 



00 

+ 



pk y (p + kY {p + kY + 6^ j 
Ret++(2+;0,p)Ret+°(z+;fc,fco(2+)) -Ret++Ret+°] (B6) 
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Finally, the matrix element containing many-Coulomb photon exchanges, is given by the 
following expression 



Im (^ol {t{Eo) g"-^(^o) t{Eo)y^ l^c 



M^okoiEo) 2 



+ l++0|2 



4:71 



(B7) 



Adding things together, and evaluating the decay width F = — 2Imz, we arrive at 
Eq. ( [4 .131) where the quantity A is given by Eq. ( [4.26|) . 



APPENDIX C: MATCHING OF THE POTENTIAL IN THE TWO-CHANNEL 

SCATTERING CASE 

In order to obtain the explicit form of the functional J^[to] in the matching condition for 
the potential (|4.36| ), we first evaluate the integrals containing go — go, which appear in the 
r.h.s. of Eq. { WM) 



Re{to(go - go)to}"^° 



thr 47r2 

Re {to(go - go)io(go - go)io}+o 



K{g^[to]-2Ml+Q^[t^]] + o{5) 



thr IGtT^ 

where t^^ = '^(2;+; 0, 0) the quantity k is defined after formula ( [4.17D , and 



(CI) 



g,[to] = J rfpt+"(2;+;0,p)C(^+;P,0), 





(C2) 



Substituting these expressions into the matching condition for the amplitude in the 
charge-exchange channel, we finally obtain 



47r2 

in , , tg 



167r2 



IT 



Q[k 



47r "'^ ' 47r3 

where the functionals Q[to] and '^[6;to] are defined in Appendix |B]. 



(C3) 
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FIGURE CAPTIONS 

FIG. 1^. Representative diagrams that contribute to the non-relativistic scattering ampli- 
tude in the purely strong case: a) tree diagram with the non-derivative coupling go; b) tree 
diagram with the derivative coupling in the four-pion vertex {gi or (72), and with mass in- 
sertions; c) one strong bubble with non-derivative couplings and with mass insertions in the 
external and internal lines. The filled circle and filled square denote the non- derivative and 
derivative vertices, respectively, and crosses stand for the mass insertions. 

FIG. |2|. Building blocks for the non-relativistic scattering amplitude at 0{a): a) one 
Coulomb photon exchange; b) strong 4-particle vertex; c),d) exchange of Coulomb photon 
in the initial and final states; e) two-loop diagram with the exchange of Coulomb photon in 
the intermediate state. 

FIG. ^. Decomposition of the relativistic scattering amplitude into one-photon exchange 
part (a) and the truncated amplitude (b). 

FIG. 1^. Summation of the mass insertions in the non-relativistic pion propagator. 

FIG. Direct quark-photon interactions in QCD and in ChPT. Solid and dashed lines 
denote vr^ and 7r°, respectively. 

FIG. ^. A particular diagram corresonding to the virtual photon correction to the scattering 
amplitude vr+vr" — > vr^vr^. Solid and dashed lines denote n"^ and vr", respectively. Filled dot 
stands for the non-minimal photon coupling to four pions. 
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